A Rigorous Proof of Fermi Liquid Behavior 
for Jellium Two-Dimensional Interacting Fermions 
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Using the method of continuous constructive renormalization group around the Fermi surface, it is 
proved that a jellium two-dimensional interacting system of Fermions at low temperature T remains 
analytic in the coupling constant A for |A|| logT| < K where K is some numerical constant and T is 
the temperature. Furthermore in that range of parameters, the first and second derivatives of the 
self-energy remain bounded, a behavior which is that of Fermi liquids and in particular excludes 
Luttinger liquid behavior. Our results prove also that in dimension two any transition temperature 
must be non-perturbative in the coupling constant, a result expected on physical grounds. The 
proof exploits the specific momentum conservation rules in two dimensions. 
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Conducting electrons in a metal at low temperature are 
well described by Fermi liquid theory (FLT). However, it 
is known that the FLT is not stable at strict zero tem- 
perature, even if the interaction is repulsive |Q. Indeed, 
below some critical temperature the effective electrons 
or holes which are the excitations of the Fermi system 
bound into Cooper pairs and the metal becomes super- 
conducting. 

With the discovery of the high temperature supercon- 
ductors (HTSC), the Fermi liquid paradigm underlying 
the BCS theory has been questioned ||]. In fact, experi- 
ments showed that the physical properties of the normal 
phase of the HTSCs are in conflict with FLT (for a recent 
experimental survey, see Ref. ||). For this reason, new 
theoretical proposals have been done by different groups 
j|. In particular, a deeper understanding of FLT itself 
became necessary. 

During the last ten years a program has been de- 
signed to investigate rigorously Fermionic many-particle 
systems by means of field theoretical methods . In 
particular the Wilson renormalization group has been ex- 
tended to models with surface singularities such as the 
Fermi surface. The ultimate goal is to create a math- 
ematically rigorous theory of the BCS transition and 
of similar phenomena in solid state physics. This is 
a long and difficult program which requires to glue to- 
gether several ingredients, in particular renormalization 
group around the Fermi surface and spontaneous sym- 
metry breaking. 

In this paper we give a contribution towards a rigor- 
ous scenario for the description of the normal phase in 
Fermionic many-particle systems. As we have already 
mentioned in the opening paragraph, FLT is not valid 
at zero temperature due to the superconducting instabil- 
ity. Even when the dominant electron interaction is re- 
pulsive, the Kohn-Luttinger instabilities JlJ prevent the 
FLT to be generically valid down to zero temperature. 
There are nevertheless two paths towards a mathemati- 
cally rigorous FLT. The superconducting instabilities can 
be avoided by considering models in which the Fermi sur- 



face is not invariant under p — ► — p Q . In two dimensions 
it is possible to prove (even non-perturbatively) that in 
this case the FLT scenario remains valid at zero temper- 
ature. However, this analysis requires a rigorous control 
of the stability of a non-spherical Fermi surface under 
the renormalization group flow, a difficult technical issue 
|Ufl| . The other alternative is to study the FLT at finite 
temperature above the superconducting transition tem- 
perature. A system of weakly interacting fermions has 
an obviously stable thermodynamic limit at high enough 
temperature, since the temperature acts as an infrared 
cutoff on the propagator in the field theory description 
of the model. In this point of view, advocated in Ref. 
|llf , the non-trivial theorem consists in showing that sta- 
bility (i.e. summability of perturbation theory) holds 
for all temperatures T satisfying to |A|| logT| < K, and 
that in that range of temperatures and couplings the first 
and second derivatives of the self-energy with respect to 
the momenta remain uniformly bounded. This is what 
we have recently completed. We developed a rigorous 
continuous renormalization group method Jl2| , and ap- 
plied it to the simplest two dimensional model of inter- 
acting Fermions, the jellium model |ll|. (This model 
is a good physical approximation of the Hubbard model 
in the regime of weak filling). In this way we proved 
a mathematically rigorous theorem, stated below, which 
physically means that the system above the BCS temper- 
ature is a Fermi liquid. Here we only outline theproof of 
this theorem, and for more details we refer to p3| . 

Let C(k) be the Fermionic propagator of the isotropic 
jellium model at finite temperature T: 
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where a, 6 € {1, 2} are the spin indices and T](k) is some 
ultraviolet cutoff (the infrared cutoff is given by the tem- 
perature). The vector k is two-dimensional. The pa- 
rameters m and /i correspond to the effective mass and 
to the renormalized chemical potential (which fixes the 
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Fermi energy), and the Matsubara frequencies kg take 
the discrete values 



2n+ 1 



n G 



(2) 



where (3 = 1/T, and we fixed 2m — 1. 

Let the system of interacting Fermions be defined by 
the Grassmann measure: 



(3) 



where cfyxc(V'jV') is the Grassmann Gaussian measure 
with covariance C given by Eq. (Q) and the interaction 
in a finite volume V is 
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(4) 

A being the bare coupling constant and Sfi being the bare 
chemical potential counterterm. The self-energy E(/c) is 
defined as the sum over all non-trivial 1PI two point sub- 
graphs (i.e those not reduced to a sin gle Su counterterm). 
The following theorem, advocated in is proved in |p_3|| : 

Theorem. The vertex functions (and the connected cor- 
relation functions) of the measure (Q) are analytic in the 
bare coupling constant X, for \X\ < c, with c given by the 
equivalent relations 
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(5) 



for some constant K (this relation is limited to the in- 
teresting low temperature regime T < 1). Furthermore in 
that range all the first partial derivatives of the self energy 
on the Fermi surface and all the second partial derivatives 
of the self energy anywhere in momentum space (in every 
direction) remain bounded: 



dki |fco = £,e(fc)=o 
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where i and j take values 0,1,2, and K\, Ki are some 
constants. 

This theorem rules out Luttinger liquid behavior, since 
for a Luttinger liquid the first derivative of E would be 
bounded only by -?GJA| in the range considered, and the 
second derivative of E would behave as A 2 /T, which is 
certainly not uniformly bounded for |A||logT| < K. In 
fact in dimension d = 1 this Luttinger liquid behavior 
has been proved up to T — |l^,|l5|. Remark that for 
the "Fermi surface" conditions we took ko 







since we 



cannot take kg = which would violate antiperiodicity 
of the Fermions. 

In the following we will give an idea of the mathemat- 
ical techniques used to obtain this result, in the hope 
that they might be of interest to a general audience of 
physicists. 

Fermionic path integrals are based on the rules of non- 
commutative or Berezin integration. This means that 
physical quantities can be expressed in terms of determi- 
nants. The correlation functions are then defined by 



S(yi, x p , z\, Zp) 



E—r — 7T / dx x ... dx n+n i det({xi},{yi},{zi}) (8) 



where dfic is the Fermionic Grassmann Gaussian 
measure with propagator C and the determinant 
det({xi}, {j/i}, {zi}) corresponds to the sum over all Wick 
contractions of the fields. 

To compute intensive quantities in the thermodynamic 
limit, one needs to sum over connected graphs only. Now, 
the simplest structures connecting points are trees, and 
each connected Feynman graph contains a tree. So why 
not reorganize the series by summing over trees only, 
keeping the remaining loop structure of the graph into 
the closed form of a determinant? This is the point of 
view adopted in Ref. JfJ]]. It has the advantage that 
since there are much less trees than graphs, and since 
determinants, by Gram's or Hadamard inequalities can 
be easily bounded, the convergence of the Fermionic se- 
ries follows. Two difficulties have to be solved. The first 
problem comes from the fact that there are typically sev- 
eral trees in a Feynman graph, so which one should one 
pick? This is a combinatoric problem which is completely 
solved by the invention of "tree formulas" which express 
in the most symmetric way how picking a tree affects the 
combinatoric of the remaining loop lines p6[ . The key 
point is that these tree formulas preserve positivity, hence 
Gram's or Hadamard inequalities 17 ll|. The connected 
amputated Schwinger functions, in the thermodynamic 
limit, may then be written as 
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where 4>i, ■■•<f>2p are some test functions we have inserted 
instead of the amputated external propagators, T n+n > is 
the set of trees connecting the n + n' vertices (hence with 
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n + n' — 1 lines) , the parameter Wi is associated by the 
tree formula to the tree line propagator Ci i , deti oop ({w i }) 
denotes the determinant of the matrix containing the un- 
contracted loop fields, whose propagators depend from 
the u>i parameters too. Without performing renormal- 
ization one can prove, for high enough temperature, that 
this sum is absolutely convergent. The strategy of the 
proof is the following. 

The determinant is bounded using a Gram inequality 
which states that 

<!<•' /<•//.. > i < n ii^n n 11*11 ( io ) 

* 3 

where fi and gj are any sets of n vectors in an Hilbcrt 
space. The spatial decay of the tree line propagators 
is used to perform the spatial integrals over X2---x n + n i , 
the first integral being performed with the help of the 
test function fa; the result of these operations is some 
function /(T, T). The sum over trees is simply bounded 
by 

E/( T ) <* T n+n> SU P /(T). (11) 
T T 

Now, by the Cayley's theorem, the number of trees sat- 
isfies #T„ +n / < (n + n') n+n and is compensated by the 
factors 1 / n\n'\ in Eq. (||) . We obtain that the series Eq. 
is bounded by 

\T2 P (fa, -fa P )\ < Y.(\AKi) n {\5ii\K 2 Y' , (12) 

where K\ and K% depend on T. This sum is convergent 
for A and 5fi small enough, depending on T. 

The tree formalism then allows us to compute the ther- 
modynamic limit of Fermionic systems with cutoffs easily 
in a transparent and convergent way. It is the right com- 
promise between the unexpanded Berezin determinant 
which is good for convergence of the series in a finite 
volume but unadapted to thermodynamic limit, and the 
Feynman graphs which are good for thermodynamic limit 
but too many for convergence of the series. 

The second difficulty was raised in particular in Ref. 
pT| : is it possible to perform renormalization and com- 
pute continuous renormalization group flows in such a 
tree formalism? This is not obvious because it is known 
that renormalization is due to the loop structure of Feyn- 
man graphs, so how to see this loop structure if loop lines 
remains unexpanded in a determinant? The answer to 
that question is twofold. First, the w parameters in the 
tree formula can be combined with a scale analysis to 
ensure that the tree which is picked in a graph is always 
optimal in the sense that hard "short-range" propagators 
are always selected in priority with respect to softer ones. 
As a result, each Wi parameter fixes the energy of the tree 
line li and constrains loop lines connecting T and T[ to 
be softer than li (where T and T[ are the two subtrees 
connected by li). Then the energy parameters Wi can be 



ordered, cutting the impulsion space in a set of bands. 
Second, once this scale analysis has been performed, a 
partial expansion of the loop determinant can detect all 
the dangerous two and four point functions which require 
renormalization. A key point is that this expansion can 
be done without destroying the Gram bound, and the 
corresponding sum is not too big (this means its cardinal 
remains bounded by K n (where K is a, constant)) be- 
cause in typical graphs there are not many two and four 
point subgraphs. 

Let us see in more detail the technique. The renor- 
malization group analysis of the theory is performed, in 
momentum space, according to the distance from the 
Fermi surface fco = tt/A e(fc) = 0. The propagator is 
decomposed over an auxiliary parameter a, running from 
1 to oo, such that this distance is roug hly or 1 ! 2 . (This 
decomposition is written as an integral to preserve the 
continuity of the renormalization group flows) . The cor- 
responding propagator C a decays on a range y/a. But 
this naive scale analysis does not lead to the right power 
counting; it is necessary to perform a second decomposi- 
tion with respect to the angular direction of k. Combin- 
ing this angular decomposition into "sectors" with the 
momentum conservation rule in two dimensions we re- 
cover the correct power counting of a just renormalizable 
field theory . 

However some complications are worth being men- 
tioned. The consideration of almost collapsed vertices 
in which the four angular directions are almost parallel 
forces us to adopt "anisotropic sectors" which have length 
Q,- 1 / 4 m t ne tangential direction, much longer than their 
radial width or 1 / 2 . 

Since a vertex is hooked to four different propagators 
with possibly different scales, the angular decomposition 
cannot be done just once, but must be updated progres- 
sively, taking also into account the momentum conserva- 
tion of larger and larger subgraphs as the renormalization 
group flow progresses towards the Fermi surface. This is 
one of the most delicate technical points in Jl3| . 

Once the correct power counting of a renormalizable 
theory has been obtained, it remains to perform correctly 
the renormalization of the Fermi radius induced by the 
interaction; indeed since the two point subgraphs are rel- 
evant rather than marginal, our bound for the critical 
temperature remains a power of A instead than exponen- 
tially small, as long as this renormalization is not per- 
formed. 

As we have already mentioned, it is possible through 
some auxiliary expansion of the loop determinant to test 
for the presence or absence of two point subgraphs to 
be renormalized, since this auxiliary expansion does not 
generate any factorial p| . But a last difficulty has to be 
overcome. Constructive renormalization has to be per- 
formed by subtracting in direct space, not in momen- 
tum space. Let G(x, y) be the amplitude of a two point 
subgraph to renormalize, and C a {. 1 x)C a {y, .) be its two 
external propagators, in an angular sector a. In the sub- 
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traction process, the tree decay inside the two point sub- 
graph must be used to compensate for the gradient term 
V — a acting on the external propagators of the subgraph. 
The tree decay inside G is anisotropic, and can happen 
in a different direction than the sector a. In that case a 
further auxiliary expansion has to be performed to take 
into account the loop structure of G(x,y); without that 
auxiliary expansion, renormalization would not work as 
expected. This last difficulty was never recognized be- 
fore. 

The bound on the first derivatives of the self-energy 
(^J) is then an easy consequence of this additional loop 
analysis. The factor A 2 in this bound already excludes 
the Luttinger liquid behavior. Finally to obtain that the 
second derivatives of the self-energy E remain bounded, 
hence to prove the Fermi liquid behavior, we have to per- 
form a second auxiliary analysis which refines anisotropic 
sectors of one four point function into isotropic sectors. 
This bound consumes a second logarithm because of the 
possibility of collapse for this four point function, and 
this explains that there is no coupling constant left in 
the bound (0). 

In three dimensions such a detailed analysis with angu- 
lar decomposition of the propagator is no longer possible 
because the momentum conservation rules are less rigid 
and allows for the possibility of non-planar vertices. But 
using the Hadamard bound on determinants (which is 
worse than the Gram bound in the sense that it con- 
sumes the factor 1/nl of symmetry), we hope to prove 
the same exponentially small bound on the BCS tran- 
sition temperature. Essentially one should combine the 
result of [[19| with a multiscale cluster expansion, and the 
renormalization of two-point functions. 

The theorem in ||l3|| , combined with four point function 
renormalization and the infrared analysis of the Gold- 
stone boson, may also lead in the future to a rigorous 
analysis of the superconducting phase of two-dimensional 
Fermi liquids. The analysis performed in []l3f should also 
be generalized in the future to study other types of insta- 
bilities. Particularly interesting is the case of a nesting 
of the Fermi surface in the Hubbard model at half-filling. 
Nesting in this case is responsible for the appearance of 
antifcrromagnetic instabilities. A mathematically rigor- 
ous analysis of this case is however very difficult and is 
not simply a trivial extension of the analysis performed 
in this paper. The point is that the diamond shaped 
Fermi surface introduces new technical difficulties as, for 
instance, the presence of van Hove singularities and flat 
sides of the diamond M| |pof . 

We thank M. Salmhofer for explaining us the mean- 
ing and importance of bounds (^0) on the self-energy 
derivatives in order to distinguish between Luttinger and 
Fermi liquid behavior. We also thank F. Nogueira for 
discussions on the content of this letter. 
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